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Complex coordinate scaling (CCS) is used to calculate resonance eigenvalues and eigenstates for 
a system consisting of an inverted Gaussian potential and a monochromatic driving field. Floquet 
eigenvalues and Husimi distributions of resonance eigenfunctions are calculated using two different 
versions of CCS. The number of resonance states in this system increases as the strength of the 
driving field is increased, indicating that this system might have increased stability against ionization 
when the field strength is very high. We find that the newly created resonance states are scarred on 
unstable periodic orbits of the classical motion. The behavior of these periodic orbits as the field 
strength is increased may explain why there are more resonance states at high field strengths than 
at low field strengths. Close examination of an avoided crossing between resonance states shows 
that this type of avoided crossing does not delocalize the resonance states, although it may lead to 
interesting effects at certain field strengths. 
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I. INTRODUCTION 

The study of time-periodic quantum systems has at- 
tracted considerable interest in recent years. One of the 
primary motivating factors for this interest is the devel- 
opment of ultra-high intensity lasers, which can produce 
electric fields within atoms that rival those produced by 
the atomic nucleus. Experiments with these ultra- intense 
lasers have led to the discovery of many interesting new 
phenomena, such as high harmonic generation Q]. Sim- 
ple one-dimensional models of the interaction between in- 
tense lasers and atoms have been shown to reproduce, at 
least qualitatively, many of these interesting phenomena 
. These models are especially interesting because their 
classical versions display chaotic motion H. In addition 
to providing insight into recent experiments, the study 
of these models can also provide insight into quantum- 
classical correspondence. 

One of the more interesting phenomena observed in 
these systems is the stabilization of atoms in intense laser 
fields. Stabilization is characterized by a decrease in the 
probability for an electron to ionize as the laser intensity 
is increased. This effect was first discovered in theoretical 
studies of the interaction between high-frequency lasers 
and atoms |I| , but recent experiments have verified that 
stabilization occurs in real atoms (|||. Studies of the 
underlying classical dynamics of these systems using one- 
and two-dimensional models have shown that the classi- 
cal motion can often account for the increased stability 
of the atom at higher laser intensities ||]7]] . 

The study of time-periodic quantum models is usually 
carried out within the context of Floquet theory || . Flo- 
quet eigenstates are eigenstates of the one-period time 
evolution operator and are the natural states for describ- 
ing time-periodic systems. In some cases the Floquet 
states of the system can be localized on stable structures 



in the classical phase space || , and this can lead to stabi- 
lization because these Floquet states have very long life- 
times. In this case, stabilization would also be predicted 
by the classical dynamics. However, there are often sig- 
nificant differences between the classical and quantum 
dynamics of chaotic systems. One of the most striking 
examples of this is scarring, where quantum eigenstates 
have higher probability to be found near the locations 
of unstable periodic orbits in the classical phase space 
JlOf . The scarring of Floquet states on unstable periodic 
orbits might make it possible for a quantum system to 
exhibit stabilization even when the corresponding classi- 
cal dynamics is unstable. Some earlier studies indicate 
that stabilization can be associated with states that are 
scarred on unstable or weakly stable periodic orbits jll) . 

In this paper we examine a time-periodic system with 
one space dimension which shows signs of stabilization. 
In this system the number of localized Floquet states, or 
resonance states, increases as the intensity of the driving 
field is increased. In Sec. [n]we present the model and dis- 
cuss the classical dynamics as well as prior studies of the 
quantum dynamics. In Sec. |ni] we describe two different 
versions of complex coordinate scaling and compare their 
predictions in this system. In Sec. |W| we investigate the 
relationship between the resonance states and the classi- 
cal dynamics of the system. We find that the resonance 
states that are created as the driving field is increased 
are associated with unstable periodic orbits in the clas- 
sical dynamics. We give an explanation, based on this 
association with periodic orbits, for why the number of 
resonance states increases as the laser field is increased. 

In Sec. [v| we carry out a detailed study of an avoided 
crossing between two resonance states. Avoided cross- 
ings in time-periodic quantum systems lead to significant 
changes in the structure of the Floquet states p"2| , |l3| |. 
Overlapping avoided crossings can even lead to delocal- 
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ization of Floquet states |13|| . Since stabilization depends 
upon the Floquet states remaining localized near (stable 
or unstable) classical structures, avoided crossings may 
play an important role in destroying stabilization. Fi- 
nally, in Sec. VI we summarize our findings. 



II. DRIVEN INVERTED GAUSSIAN MODEL 



The model we will study is an inverted Gaussian po- 
tential interacting with a monochromatic driving field in 
the radiation gauge. The Hamiltonian of the system in 
atomic units (which are used throughout the paper) is 



ff=~(p-~sin(wt)) -i; )( . xl )(-( ,-.„)-'] (I) 



where Vo = 0.63 a.u. and a = 2.65 a.u. It is useful to 
write this as H = Hq + V, where 



H 



Vbexp(-(a;/a) 2 ) 



and 



V 



-p sin(wi) 



2uj- 



■ sin 2 (u;£). 



(2) 



(3) 



Figure [l] illustrates the classical dynamics of this system 
for driving frequency ui = .0925 a.u. The strobe plots in 
Fig. [j] are calculated by evolving a set of trajectories, all 
with initial momentum p — 0, over many cycles of the 
field and plotting the location of each trajectory when 
t = Inn j io (after each full cycle of the field). When no 
driving field is present the motion is regular and bounded 
for negative energies. Motion at positive energies is un- 
bounded. Figs. [l]a, [l]b, and 0c show the classical strobe 
plots for e = .038, .065, and .09 a.u., respectively. As e 
is increased the region near (x = 0,p = 0) remains sta- 
ble, but the size of the stable region gets smaller as e is 
increased. The filled squares in Fig. [I] indicate the lo- 
cations of the periodic points in the strobe plot and the 
arrows show x = a and x = 2a, where a = e/u> 2 is the 
excursion parameter of a free electron in the field. At 
these parameter values the periodic orbit at (0, 0) is sta- 
ble while the other two periodic orbits are unstable. As 
e is increased the unstable periodic orbits move toward 
larger values of x. For very high frequency driving fields 
two of the periodic orbits can be stable while the third is 
unstable. This is illustrated in Fig. [l]d, which shows the 
classical strobe plot for to = 2 a.u. and e = 42 a.u. The 
value of a in Fig. |l|d is the same as in Fig. [ljs, but at the 
higher frequency the periodic orbit located near x = 2a 
is a stable elliptic orbit surrounded by regular motion. 
The periodic orbit at x — a is hyperbolic. 
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FIG. 1. Strobe plots of the classical dynamics for the driven 
inverted Gaussian system. The initial conditions used to gen- 
erate the plots all lie on the line p — 0. a is the classical 
excursion parameter for a free electron the field. The loca- 
tions of the periodic orbits (stable and unstable) are indicated 
by filled squares. 



The quantum dynamics of this system has been the 
subject of several investigations during the past decade. 
The resonance states of this system were first calculated 
by Bardsley and Cornelia in 1989 (m). More recent stud- 
ies have focused on high-harmonic generation (HHG) in 
this system It is the findings of Ben-Tal, Moi- 

seyev, and Kosloff jl^j , hereafter BMK, that have the 
most relevance to our work. They found that the num- 
ber of resonance states in this system increased as the 
field strength was increased over a certain range. BMK 
attempt to explain the creation of new resonance states 
as the field strength is increased by analyzing the dy- 
namics of the time-averaged system in a reference frame 
that oscillates with a free electron in the driving field, 
known as the Kramers-Henneberger or K-H frame @- 
They found qualitative agreement in that the number of 
bound states in the time-averaged potential increases as 
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the field strength is increased. However, the quantita- 
tive agreement was not very good. This is not surpris- 
ing since the time-averaged K-H description is only accu- 
rate for very high frequency driving fields. Since the fre- 
quency of the driving field used by BMK and in this work 
(u> = .0925 a.u.) is lower than the frequency of motion 
for two of the bound states in the undriven system (.4451 
a.u. and .1400 a.u.), the time-averaged K-H description 
is not quantitatively accurate. It is somewhat surprising 
that the time-averaged K-H description is qualitatively 
accurate because the classical motion of the system in 
the time-averaged K-H frame is stable while the classi- 
cal motion of the true system is largely unstable. For 
a > 1 the time-averaged potential in the K-H frame is 
a double well with minima separated by approximately 
2a. Motion in this double well would be quite different 
from that seen in the strobe plots in Fig. [j]a-c (although 
it would closely resemble the motion shown in Fig. [yd, 
which is at a frequency that is high enough for the time- 
averaged K-H description to be valid). Our goal in this 
paper is to find an alternative explanation for the cre- 
ation of resonance states as e is increased in this system, 
an explanation that does not rely on the time-averaged 
K-H description. 



III. COMPLEX COORDINATE SCALING 



integrable eigenfunctions. As a result of this scaling the 
new time-independent Hamiltonian is 



H = H (xe ie ) = tl^L - ^ exp (-(xe^/a)) . (5) 



The kinetic energy operator is easily evaluated using the 
basis states in Eq. B As l ong as our box is sufficiently 
large (L » 2xq/ y/cos(28)) we find that 



(m\ - Vo exp (-{xe l9 /a) 2 ) \n) = V(m + n) - V(\m - n\) 

(6) 



where 



V(j) = exp 



U? ) C ° S ( "2~ 



(7) 



In recent years the technique of complex coordinate 
scaling (CCS) has been used extensively in the study of 
open quantum systems. In this section we will review 
two versions of complex coordinate scaling (standard and 
exterior scaling) and show how these techniques can be 
used to compute the resonance states of an open, time- 
periodic system. Results from the standard and exterior 
scaling versions are compared, for both time- independent 
and time-dependent calculations. 



A. Standard complex coordinate scaling (CCS) 



Once these matrix elements are calculated the Hq matrix 
can be constructed. Diagonalizing Ho yields the energy 
eigenvalues of the time-independent system as well as the 
eigenvectors 



N 

£■ 

71=1 



(8) 



We first examine how the eigenvalues and eigenstates 
of a time-independent open system can be calculated us- 
ing standard complex coordinate scaling (CCS), a tech- 
nique that is examined in detail in |ilj|Jl9| |. In this paper 
we will use a basis of particle-in-a-box states for our cal- 
culations. These states arc defined by 

J~2 . (u-kx nir\ 

where —L/2 < x < L/2. Calculations using CCS are per- 
formed just as they are in traditional quantum mechan- 
ics, except that the coordinate is scaled in the Hamilto- 
nian so that x — > xe lS (0 < 6 < tt/4). Scaling the co- 
ordinate in this fashion allows us to represent resonance 
states, which are not in the Hilbert space, using square 



Fig. £p shows the energy eigenvalues of Hq calculated 
without complex scaling (8 = 0). The potential sup- 
ports three bound states at E = —0.4451, —0.1400, and 
—0.00014 a.u. Without complex scaling all eigenvalues lie 
on the real axis. When the coordinate is scaled, though, 
the Hamiltonian becomes non-Hermitian and it is pos- 
sible for eigenstates of the scaled system to have com- 
plex eigenvalues. This can be seen in Fig. ||b, which 
shows the eigenvalues calculated using CCS with = 0.3. 
The bound state eigenvalues remain on the real axis but 
the positive-energy continuum states are rotated into the 
lower half plane by an angle of 28. It is this rotation of 
the continuum that will allow us to identify resonances. 
No resonances exist for the system with Hamiltonian Hq , 
only bound states and a continuum of energy eigenvalues. 
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FIG. 2. The complex-scaled energies of the undriven in- 
verted Gaussian system. The unsealed energies are shown in 
(a), the CCS energies in (b), and the ECCS energies in (c). 
The bound states of the system have energies -0.4451, -0.1400, 
and -0.0001 a.u. All calculations were performed using a box 
size L — 200 a.u. and 400 basis states. 



B. Exterior complex coordinate scaling 

The basic idea of exterior complex coordinate scaling 
(ECCS) is to scale the x coordinate by a factor e l9 as in 
CCS, but only in the region \x\ > x s where the potential 
is zero. Discontinuities at ±x s are avoided by using a 
smooth scaling relation x — ► F(x), where 



F(x) = x + {e ld - 1) 



1 / cosh(A(x — x a )) 
2A \cosh(A(x + x s )) 



(9) 



with A = 5 a.u. and x s — 25 a.u. This exterior scal- 
ing method is given a thorough presentation in [ ]l9|j20|] . 
Because the potential is zero in the region where the co- 
ordinate is scaled, the potential matrix elements can be 
calculated without any complex scaling (i.e. using Eq. 
^ and [?], but with = 0). The scaled time-independent 
Hamiltonian then becomes Hq = Hq + Vcap , where 



V CA p(x) = Vo(x) +Vi(x)^ 



d 2 



(10) 



acts as a complex absorbing potential. The coordinate- 
dependent factors in Eq. \L0\ are defined by 



V (x) 



0£ 

dx 



(11) 



(12) 



and 



v 2 {x) = -(i-r 2 {x)), 



(13) 



where f(x) = dF/dx. Plots of F(x), V (x), Vi(x), and 
V^(x) are shown in Q]2(| . 
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FIG. 3. Wavefunctions of two ECCS continuum states. 
The state shown in (a) is a partially-scaled state whose eigen- 
value is rotated by less than 29 from the real axis. The state in 
(b) is a fully-scaled state whose eigenvalue is rotated the full 
28 from the real axis. The partially- scaled state is localized 
between — x 3 and x s (x a = 25 a.u.), while the fully-scaled 
state is almost excluded from this region. Both states are 
peaked near ±x a . 

We will again use a basis of particle-in-a-box states to 
calculate Hq. Matrix elements for the potential energy 
term are calculated without any complex scaling, while 
the kinetic energy and Vcap matrix elements are calcu- 
lated numerically. Diagonalizing Hq gives the complex 
energy eigenvalues for the exterior scaled system, which 
are shown in Fig. ^c. Note that the bound state eigen- 
values are still on the real axis and most of the contin- 
uum states have been rotated into the lower half plane by 
26. However, several of the positive energy states have 
been rotated into the lower half plane by considerably 
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less than 29. We refer to these as "partially scaled" con- 
tinuum states. Fig. [| shows the wavefunction of one fully 
scaled continuum state and one partially scaled contin- 
uum state. The partially scaled state is strongly peaked 
near x = x s and it is non-zero only within the region 
—x s <x<x s , while the fully scaled state is zero within 
this region. As x s is decreased toward 0, the number of 
partially scaled states decreases. At x s = the ECCS 
eigenvalues exactly match the CCS eigenvalues as ex- 
pected. 

Since we will want to examine the structure of the res- 
onance wavefunctions in the periodically driven system it 
is important first to examine the structure of the eigen- 
states of Hq. We will examine the structure of the eigen- 
states of Hq by calculating Husimi distributions for 
each of the three bound states. An Husimi distribution is 
a quasi-probability distribution of a quantum state in the 
phase space. The Husimi distribution (HD) of a quantum 
wavefunction Hf(x) is defined as 
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FIG. 4. Husimi distributions of the bound states of the in- 
verted Gaussian system. The distributions in the left column 
are of the CCS bound states. The distributions in the center 
column are CCS bound states rotated back into the real co- 
ordinate frame. The right column shows the distributions for 
the ECCS (x s — 25 a.u.) bound states. The scaling angle is 
9 = 0.3 for all plots. 

Husimi distributions for the three bound states of Hq 
are shown in Fig. || The left-hand column shows the 
HDs of the eigenstates, ipi(x), of Hq obtained using CCS. 
The center column shows the HDs for the eigenstates 



i\)i{xe %e \ which is ipi(x) rotated back to the real coordi- 
nate frame [p2| . The right-hand column shows the HDs 
of the eigenstates that were calculated using ECCS. Note 
that the center and right-hand columns agree quite well, 
while the left-hand column shows significant differences. 
The Husimi distributions in the right-hand column match 
those that are found without complex scaling the Hamil- 
tonian and hence they are the correct distributions. The 
distributions in the center column are correct except for 
the patches of black in the top-left and bottom-right cor- 
ners. These patches show where the Husimi distribution 
blows up because the wavefunction is not square inte- 
grable. 



C. Floquet calculations 

As we have seen, complex coordinate scaling can be 
used to calculate the first N energy eigenstates of Ho. 
These eigenstates can then be used as a basis to com- 
pute the one-period time evolution (Floquet) matrix, 
C/(T), for the driven system. This matrix is calculated by 
numerically integrating the time-dependent Schrodinger 
equation TV times from t = 0tot = T = 2n/uj with initial 
conditions \^f(t = 0)) = l^i), where is the ith energy 
eigenstate of Hq. Diagonalization of this matrix gives 
the Floquet eigenvalues and eigenstates (in the basis of 
eigenstates of Ho). 

The time-dependent Schrodinger equation for the 
driven inverted Gaussian system is 

d ~ e e 2 

ih—m = H \V) - ~psin(wt)|¥} + -^sin 2 (^)|*) 
at uj 2u> z 

(15) 

where p is the complex scaled momentum operator. Since 
all computations are performed in a basis of eigenstates 
of Hq we must first calculate the matrix elements of p. 
In calculating these matrix elements it is critical to rec- 
ognize that Ho is not a Hermitian matrix and thus its 
eigenvectors do not have the usual properties that eigen- 
vectors of Hermitian matrices have. One cannot obtain 
the left eigenvectors of a non-Hermitian matrix simply by 
taking the complex conjugate of the right eigenvectors. 
In our case, Hq is complex symmetric and the coefficients 
of the left eigenvectors are equal to (not complex conju- 
gates of) the coefficients of the right eigenvectors, so 



(16) 



n=l 



The normalization of the eigenvectors is also differ- 
ent. For our complex symmetric matrix the eigenvectors 
should be normalized so that the sum of the squares of 
the c n j's is 1, rather than the sum of the absolute squares. 
With this in mind we can calculate the matrix elements 
for p using 
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{i>i\p\ipj 



N N 

EE< 

m— 1 n— 1 



(m|p|rc) 



(17) 



The (m\p\n) are easy to calculate when CCS is used and 
p is simply pe~ %e . However, when ECCS is used those 
matrix elements are calculated numerically using 



(m\p\n) = 



ihim 



[P(m + n) + P(m - n)] (18) 



where 



P(k) = l^sin^-^) f-\x)dx 



(19) 



and f(x) is defined in Sec. [II B. 

By numerically integrating the Schrodinger equation 
over one cycle of the driving field we can construct the 
Floquet matrix as described above. This matrix is then 
numerically diagonalized to give the Floquet eigenvalues 
and eigenstates. Since the Floquet eigenstates are calcu- 
lated in a basis of eigenstates of Hq we can write them 
as 



JV 



(20) 



i=l 



Because they are eigenstates of the one-period time evo- 
lution operator (Floquet matrix) we can write 



U{T)\qp) 



"in-' 



(21) 



where qp is the quasienergy of the state {qp}- Because the 
Hamiltonian Hq is not Hermitian, the time evolution op- 
erator is not unitary. This means that the Floquet eigen- 
values do not necessarily have unit modulus, and thus the 
quasienergies qp are in general complex. We can write the 
quasienergies as qp — flp + iTp/2, where rg = 1/r^ is 
the lifetime of the state \qp). Resonance states are easily 
identified by plotting the Floquet eigenvalues, which we 
will denote as = exp(-iqpT). Fig. || shows the Flo- 
quet eigenvalues calculated using both CCS and ECCS 
for the driven Gaussian system with u> = 0.0925 a.u. and 
e = 0.038 a.u. The eigenvalues that were found with the 
CCS method form a well-defined spiral from the origin 
out to the edge of the unit circle. These states are indi- 
cated by filled circles in Fig. [|a. Resonance states are 
indicated by filled squares and lie off of the continuum 
spiral. The continuum spiral is not as well-defined when 
the ECCS method is used, as shown in Fig. [pD. How- 
ever, only a few eigenvalues near the origin appear to 
fall out of the spiral. This could cause some difficulty 
in identifying broad (short-lived) resonances, but narrow 
(long-lived) resonances can still be easily identified. Fig. 

shows that CCS and ECCS appear to give the same 
resonance eigenvalues. 
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FIG. 5. Floquet eigenvalues for the driven inverted Gaus- 
sian system with ui — 0.0925 a.u. and e = 0.038 a.u. The 
eigenvalues calculated using CCS are shown in (a). In (b) the 
ECCS eigenvalues are shown. The scaling angle is 9 — 0.3 for 
both sets of eigenvalues. Resonance states are indicated by 
filled squares, while continuum states are indicated by filled 
circles. 



The resonance eigenvalues should be independent of 
the scaling angle 9, while the continuum eigenvalues ro- 
tate around the origin as 9 is changed. However, when 
calculations are performed using a finite basis the reso- 
nance eigenvalues will be weakly dependent upon 9 [ p3| . 
To accurately determine the quasienergies (and hence the 
lifetimes) of these states it is important to optimize 9 
by finding the stationary point of each resonance eigen- 
value as 9 is changed. However, since our goal is not 
an accurate quantitative determination of eigenvalues or 
lifetimes but a qualitative understanding of the relation- 
ship between the quantum dynamics and the classical 
motion, it is not critical that 9 be optimized for our cal- 
culations. Optimizing 9 presents a problem in this type 
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of study because the optimal value of 9 is generally dif- 
ferent for different resonance states. We wish to study all 
of the resonance states of the system at the same time 
and it is impossible to optimize 9 for all resonance states 
within a single calculation of the Floquet matrix. We find 
that changing 9 between 0.3 and 0.7 results in no visible 
change in the plots of the resonance eigenvalues. There 
is also no visible change in the Husimi distribution of the 
states. Not optimizing 9 may lead to slight inaccuracies 
in the calculated lifetimes for the resonance states, but 
we find that the error in the lifetimes is no greater than 
±0.1T which is acceptable for our purpose here. 
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FIG. 6. Husimi distributions of the resonance states at 
e = 0.038 a.u. The left column shows the Husimi distribu- 
tions of the states computed using CCS. The center column 
shows the distributions of the CCS states rotated back into 
the real coordinate frame. The right column shows the distri- 
butions of the ECCS states (x s = 25 a.u.). In the right-hand 
column the locations of the periodic orbits are indicated by 
filled circles and the contour plots have been lightened to al- 
low the circles to show. The scaling angle is 9 — 0.3 for all 
plots. The lifetimes for each state are given in units of the 
driving period T = 2-k/cj. 

In Figure ^ we show the Husimi distributions of the 
three resonance states indicated in Fig. ||. As in Fig. |] 
the left-hand column shows CCS states, the middle col- 
umn shows the rotated CCS states, and the right-hand 
column shows ECCS states. Again we find agreement 
between the middle and right-hand columns, while the 
left-hand column is inaccurate (as expected). Lifetimes 
of the three states are indicated in units of the driving 
period T — 2ir/ui. Filled circles indicate the locations of 
the classical periodic orbits. The resonance state with 
the longest lifetime is almost indistinguishable from the 
ground state of the undriven system shown in Fig. ^a. 
The state shown in Fig. &> has a much shorter lifetime 



and is beginning to elongate toward the positions of the 
unstable periodic orbits, with a peak of probability near 
the periodic orbit at (x = 6.98 a.u., p = 0). The state 
shown in Fig. |^c has a very short lifetime and its Husimi 
distribution is similar to that of a continuum state. 




Re A 



FIG. 7. Floquet eigenvalues for lu — 0.0925 a.u. and two 
different field strengths. The eigenvalues are calculated using 
ECCS with 6 = 0.3 and x s = 25 a.u. Resonance states are 
indicated by filled squares and continuum states by filled cir- 
cles. At these higher field strengths the number of resonance 
states is greater than at e = .038 a.u. 



IV. RESONANCE CREATION AND SCARRING 

Because the rotated CCS states blow up for large 
|x| we choose to perform the remainder of our calcula- 
tions using ECCS only. As long as our unsealed region 
(— x s < x < x s ) is large enough to contain the phase 
space region in which we are interested, we will not have 
to worry about any inaccuracy associated with the scal- 
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ing procedures. We choose x s = 25 a.u. which is large 
enough to include the locations of the periodic orbits at 
all of the field strengths we study. 

Figure [?] shows the Floquet eigenvalues for e = 0.065 
and 0.09 a.u. Again, resonance states are indicated by 
filled squares while continuum states are indicated by 
filled circles. Comparing these plots with Fig. |^ we see 
that the number of resonance states increases as e is in- 
creased, from only three at e = 0.038 a.u. to five at 
e = 0.09 a.u. This is in agreement with BMK In 
the classical system, however, the stable structure near 
(x = 0,p = 0) gets smaller as e is increased. If the res- 
onance states were associated with this stable classical 
structure then some of the resonances should disappear 
as e is increased. Instead, the opposite behavior is found. 
To find the explanation for the increase in the number of 
resonance states we examine the Husimi distributions of 
the resonance states indicated in Fig. 
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FIG. 8. Husimi distributions for the four resonance states 
at e = 0.065 a.u. Lifetimes for each state are given in units of 
the driving period T = 2-k/lu. The locations of the periodic 
orbits are indicated by filled circles. 

Figure || shows HDs for the four resonance states at 
e = 0.065 a.u. Lifetimes for the states are indicated in 
units of the driving period and the positions of the pe- 
riodic orbits are indicated by filled circles. The state 
with the longest lifetime looks very much like the ground 
state of the undriven system. The states shown in Figs. 
||b and ||c show some similarities to the excited states of 
the undriven system, but they have both been elongated 
in the direction of the unstable periodic orbits. The state 
shown in Fig. Ijs has a probability peak near the unsta- 
ble periodic orbit at fa; = 16.67 a.u., p = 0), while the 
state shown in Fig. |8p has a peak between the two un- 
stable orbits. These two states appear to have become at 
least partially associated with the unstable periodic or- 



bits. The state shown in Figure |8jd is the newly created 
resonance and it has the shortest lifetime of the four. It 
has a modest peak near the periodic orbit at (x = 9.09 
a.u., p = 0). 
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FIG. 9. Husimi distributions for the four of the five reso- 
nance states at e = 0.09 a.u. The resonance state with the 
shortest lifetime is not shown because its Husimi distribu- 
tion is indistiguishable from that of a continuum state. Life- 
times for each state are given in units of the driving period 
T = 2n/u. The locations of the periodic orbits are indicated 
by filled circles. 

Figure |9| shows HDs for four of the five resonance states 
at e = 0.09 a.u. The state that closely resembles the un- 
driven ground state (Fig. |^b) no longer has the longest 
lifetime. Instead, the longest-lived state resembles the 
first excited bound state of Hq, but with additional peaks 
near the periodic orbits at x = 11.52 a.u. and x = 22.04 
a.u. The state shown in Fig. ^d is similar to the state 
shown in Fig. ||d, but with a more prominent peak near 
the periodic orbit at x = 11.52 a.u. Note that the life- 
time of this state is also greater than that of the state 
shown in Fig. ||d. 

At low values of e all of the resonances have their prob- 
ability concentrated near (x = 0,p — 0). At these low 
values of e the two unstable periodic orbits are located 
close to the stable orbit near (x = 0,p = 0). If any res- 
onance state was associated with the unstable periodic 
orbits at such low field strengths it would be difficult to 
tell from its Husimi distribution. As e is increased the 
unstable periodic orbits move toward larger values of x 
and some of the resonance states begin to spread in that 
direction as well. At moderate values of e some states 
show peaks near the periodic orbit that is farthest from 
(x = 0,p = 0), close to x — 2a. Only at high values of e 
do we begin to see a state that is peaked on the unstable 
orbit that is closest to (x = 0,p = 0), near x = a. We 
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believe that it is the association between the resonances 
and the unstable periodic orbits that explains the cre- 
ation of resonance states as e is increased. At low e all 
three periodic orbits are too close together to support 
many quantum states because they all occupy essentially 
the same region of phase space. As e is increased the un- 
stable periodic orbits move away from the stable one and 
from each other. This allows quantum states to be asso- 
ciated with these unstable orbits without occupying the 
same region of phase space as the states associated with 
the stable orbit, so new resonance states are created. It is 
the scarring of resonance states on unstable periodic or- 
bits of the classical system that accounts for the increase 
in the number of resonance eigenstates, even as the stable 
region in the classical phase space is diminished. 

The behavior we see here would be unlikely to sta- 
bilize the ground state of the undriven system against 
ionization in a high intensity field. This is because the 
lifetime of the resonance state that seems most closely 
related to the ground state (shown in Figs, ^a, |]a, and 
|9|b) decreases as e is increased. However, the observed 
behavior could lead to stabilization for an excited state 
of the undriven system. The excited states have most 
of their probability away from (x — 0,p — 0) and would 
thus overlap with resonance states that are not peaked at 
that point. Since these resonances grow in number and 
increase their lifetimes as e is increased, an excited state 
of the undriven system may become stabilized against 
ionization as e is increased. 



V. AVOIDED CROSSINGS BETWEEN 
RESONANCES 

Avoided crossings between resonance eigenvalues have 
been identified in this system Avoided cross- 

ings between Floquet eigenvalues play an important role 
in multi-photon ionization [^4| and the derealization 
of Floquet eigenstates In this section we investi- 

gate the quantum dynamics at one avoided crossing to 
determine if it leads to derealization of the resonance 
eigenstates. Derealization is closely related to ioniza- 
tion in these systems because long-lived resonance states 
can only exist if they are localized within the interac- 
tion region. If avoided crossings lead to derealization 
they would also lead to a decrease in the lifetime, and 
eventually the destruction, of the resonance states. 

Figure shows the Floquet eigenvalues of three reso- 
nance states at several field strengths between e = 0.076 
and 0.085 a.u. Two of these states (labeled A and B 
in Fig. and indicated by filled circles and squares, 
respectively) are involved in a prominent avoided cross- 
ing at a field strength of about e = 0.0805 a.u. The 
third resonance eigenvalue (labeled C and indicated by 
filled triangles) passes close by the other two at this field 
strength, but it is not clear from Fig. |l^ if that state is 
involved in the avoided crossing. 
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FIG. 10. Floquet eigenvalues of three resonance states (la- 
beled A, B, and C) that are involved in an avoided crossing 
around e = 0.0805 a.u. The points show the eigenvalues for 
e = .076, .078, .08, .0805, .081, .083, and .085 a.u. The num- 
bers shown in the plot indicate the field strengths at the end 
points of each eigenvalue sequence. The lines are intended 
only as an aid to the eye. 

To determine the effect of this avoided crossing on the 
resonance states we examine the Husimi distributions of 
the states A, B, and C shown in Figure As e is 
increased from 0.078 a.u. to 0.0805 a.u. the states A 
and B undergo strong mixing with each other. When 
the field strength is increased to 0.083 a.u. we find that 
states A and B have completely exchanged their struc- 
ture. State C does not appear to have any significant 
structural changes in this range of field strengths. How- 
ever, it should be noted that state C has a significant 
increase in its lifetime as e is increased from 0.078 a.u. 
to 0.083 a.u. States A and B exchange lifetimes as well 
as structure, but the lifetimes of both states at e = 0.083 
a.u. are somewhat smaller than the corresponding life- 
times at e = 0.078 a.u. It may be that state C somehow 
gains stability at the expense of states A and B, even 
though it does not appear to pick up any of the struc- 
ture of those states. Note that there are slight differences 
between the Husimi distributions of states A and B at 
e = 0.073 a.u. and the corresponding distributions at 
e = 0.083 a.u., but these may be due to a small amount 
of mixing with state C or with continuum states. 

Figure [ll] does not reveal any significant increase in 
the derealization of any of the resonance states. It does 
indicate that an avoided crossing of this type might lead 
to changes in the lifetimes of the resonance states, but 
although two states have their lifetimes decreased the 
third has its lifetime increased. There is no clear indica- 
tion that this type of avoided crossing contributes to the 
destruction of the resonance states that might prevent 
stabilization at very high intensities of the driving field. 
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We expect that the destruction of resonance states occurs 
primarily as the result of coupling between a resonance 
state and the continuum rather than between resonance 
states 



lead to increased HHG. In fact, increased HHG has been 
observed in previous studies of the avoided crossings in 



this system 
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FIG. 11. Husimi distributions for the states involved in the 
avoided crossings shown in Fig. hd. The labels A, B, and C 
correspond to the labels in Fig. jlO| . Lifetimes for each state 
are given in units of the driving period T — 2ir/uj. The top 
row (a-c) shows Husimi distributions for the three resonance 
states at e = .078 a.u. The second row (d-f) shows the dis- 
tributions at e = .0805 a.u. The bottom row (g-i) shows the 
distributions at e = .083 a.u. States A and B appear to ex- 
change their structure as they pass through the avoided cross- 
ing. State C does not appear to undergo any major changes 
in its structure, but its lifetime increases dramatically as it 
passes through the avoided crossing. The locations of the 
periodic orbits at each field strength are indicated by filled 
circles. 



Avoided crossings between resonance states may play 
an important role in other phenomena in this system. 
For example, it has been shown in bound systems that 
avoided crossings between Floquet eigenstates can result 
in increased high harmonic generation (HHG). Avoided 
crossings contribute to increased HHG in two ways. Dur- 
ing the turn-on of a laser field avoided crossings can 
put the quantum system into a superposition of Floquet 
states that may emit radiation at higher frequencies than 
would be emitted by a single Floquet state [£6[ . Avoided 
crossings also contribute to HHG by spreading the Flo- 
quet states over a wider range of energy, thus allowing a 
single Floquet state to emit higher frequency radiation. 
For the type of avoided crossing observed here the states 
are only delocalized near the exact field strength at which 
the avoided crossing occurs, because at this field strength 
the Floquet states have mixed their structure |T^]. At 
that particular field strength, though, this effect could 



VI. CONCLUSION 

We utilize the complex coordinate Floquet method 
(CCFM) to calculate resonances for an inverted Gaus- 
sian potential driven by a monochromatic field. As has 
been previously observed, we find that the number of res- 
onance states increases as the field strength is increased. 
This behavior in the quantum system seems to be op- 
posite to what is observed in the classical system, where 
the dynamics becomes increasingly unstable as the field 
strength is increased. An examination of the Husimi dis- 
tributions of the resonance states in this system shows 
that the newly created resonances states are associated 
with unstable periodic orbits in the classical motion. 
This scarring of the eigenstates on unstable periodic or- 
bits has been seen in other systems and it represents one 
of the most significant deviations of quantum dynamics 
from the corresponding classical dynamics. In this sys- 
tem, resonance eigenstates are scarred on unstable peri- 
odic orbits and the movement of the periodic orbits in 
the phase space allows for the creation of new resonance 
states as the field strength is increased. The creation of 
these new resonance states might help to stabilize the 
system against ionization in intense fields. 

An avoided crossing between resonance states is closely 
examined and it is found that, although two states do ex- 
change their structure during the avoided crossing, there 
is no significant derealization of the resonance states as 
a result of this avoided crossing. Coupling to the con- 
tinuum is more likely to destroy resonance states than 
is coupling between two resonance states. Avoided cross- 
ings of the type studied here, although they may not lead 
to faster ionization, could lead to many interesting effects 
such as increased high harmonic generation. 
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